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Abstract 

Aberrant methylation of specific CpG sites at tine promoter is widely responsible for genesis and developnnent of various 
cancer types. Even though the microarray-based methylome analyzing techniques have contributed to the elucidation of 
the methylation change at the genome-wide level, the identification of key methylation markers or top regulatory networks 
appearing common in highly incident cancers through comparison analysis is still limited. In this study, we in silico 
performed the genome-wide methylation analysis on each 10 sets of normal and cancer pairs of five tissues: breast, colon, 
liver, lung, and stomach. The methylation array covers 27,578 CpG sites, corresponding to 14,495 genes, and significantly 
hypermethylated or hypomethylated genes in the cancer were collected (FDR adjusted p-value <0.05; methylation 
difference >0.3). Analysis of the dataset confirmed the methylation of previously known methylation markers and further 
identified novel methylation markers, such as GPX2, CLDN15, and KL. Cluster analysis using the methylome dataset resulted 
in a diagram with a bipartite mode distinguishing cancer cells from normal cells regardless of tissue types. The analysis 
further revealed that breast cancer was closest with lung cancer, whereas it was farthest from colon cancer. Pathway 
analysis identified that either the "cancer" related network or the "cancer" related bio-function appeared as the highest 
confidence in all the five cancers, whereas each cancer type represents its tissue-specific gene sets. Our results contribute 
toward understanding the essential abnormal epigenetic pathways involved in carcinogenesis. Further, the novel 
methylation markers could be applied to establish markers for cancer prognosis. 
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introduction 

Carcinogenesis is a complicated cellular process starting from 
genesis of a few cancer cells and developing into malignancy 
through cellular proliferation [1,2]. Depending on the specific 
tissue and cancer stage, different sets of genes act on the cancerous 
cells in different modes [3] . For example, BRC A 1 is a well-known 
tumor suppressor gene implicated in the predisposition of early- 
onset breast and ovarian cancer [4] . On the contrary, loss of RB is 
well documented in many human tumor types [5] . 

In the course of cancer development, tumor suppressors as well 
as proto-oncogenes undergo genetic and epigenetic changes [6-9] . 
Promoter methylation is an epigenetic regulation, and its 
abnormal changes in cancer-related genes have been known to 
have pivotal roles during carcinogenesis. So far, numerous genes 
have been identified as epigenetic markers that can predict 
occurrence of specific cancers [10-13]. BRCAl is a tumor 
suppressor, and it is hypermethylated in a limited number of 
cancer types [14,15]. Meanwhile, another tumor suppressor, 
PTEN, is hypermethylated in a wide range of tumor types [16- 
1 8] . The methylation change of tumor-related genes may appear 
at a specific stage or sub-tissue. For example, methylation of 
RASSFl was observed with gradual increase in association with 



increasing tumor size and advancement of the stage of breast 
cancer [19]. In another study, ADAM33 was shown to be silenced 
by the promoter methylation with hypermethylation in invasive 
lobular carcinoma compared to invasive ductal carcinoma [20] . 

These facts imply that specific set of genes are diflFerentially 
methylated depending on the specific cancer type. Studies chasing 
methylation change have been extensively carried out in each type 
of cancer tissue. However, comprehensive approaches wherein the 
methylation profiles are compared to thus reveal the commonal- 
ities in various cancers or cancer-specific epigenetic pathways have 
been limited. Accumulation of methylation markers in cancer 
prompts us to establish a systemic approach that is able to predict 
cancer-related genes and pathways with abnormal methylation. 
Current high throughput techniques have enabled us to identify 
genes that show differential methylation in cancer tissues, and 
genome-wide methylome data are available in public databases 
such as the Gene Expression Omnibus (GEO) [21]. In this study, 
by utilizing the GEO methylome databases, we have aimed to 
elucidate the molecular pathway in which genes are highly 
dysregulated by promoter methylation and to compare them in the 
five major cancers: breast, colon, liver, lung, and stomach. To do 
this, we first analyzed in silico the lUumina methylation bead chip 
data from five cancers. After collecting data for ten pairs of normal 
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Figure 1. Cluster and heat map analysis with genes of altered methylation for five cancers. A. Cluster analysis was carried out using a 
methylation index (P-value) taken from lllumina Human Methylation 27 Bead Chip database registered at the GEO. Hierarchical clustering with the 
correlation distance was shown for 50 tumor tissues from the breast, colon, liver, lung, and stomach cancer. B. Hierarchical cluster analysis was 
conducted for 90 genes that were related with "cancer" from the five cancer types. The bottom panel is a heat map analysis wherein the rows and 
columns represent genes and cancer cell lines, respectively. 
doi:10.1371/journal.pone.0097818.g001 



and cancer tissues from each cancer type we compared the 
interrelationship of the five cancers based on methylomic change. 
Further, the Ingenuity pathway analysis was carried out to 
elucidate the top pathway that is eligible to be distorted by 
methylation. 

Materials and Methods 

Cell culture 

MCFIOA (normal breast), MCF7 (breast cancer), and HEK293 
were purchased from the American Type Culture Collection 
(ATCC; Manassas, VA). HEL 299 (normal lung) and A549 (lung 
cancer), and AGS (stomach cancer), THLE-3 (normal liver) and 



Hep3B (liver cancer), CCD-I8C0 (normal colon) and HT-29 
(colon cancer) were purchased from the Korean Culture Type 
CoUection (KCTC; Korea). MCFIOA was grown in MEBM 
supplemented with MEGM Single Qtiots and cholera toxin 
(Lonza, Basel, Switzerland). HEK293, Hep3B, THLE-3, and 
CCD-I8C0 were grown in Dulbecco's modified Eagle's medium 
(DMEM) supplemented with 10% fetal bovine serum (FBS). 
MCF7, A549, AGS, and HT-29 were cultured in RPMT1640 
supplemented with 10% FBS. 
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Figure 2. Highest confidence networl< of genes displaying altered methylation levels in breast cancer. In the network, hypermethylated 
genes In cancer are colored in red, whereas the hypomethylated genes are shown in green. The intensity of the color reflects the magnitude of 
methylation change. According to IPA, the network is relevant to "cellular movement, cellular development, and cancer." The highest confidence 
network for the other four cancers is added in Fig. S2. Each interaction is supported by at least one literature reference, with solid lines representing 
direct interactions, and dashed lines representing indirect interactions. 
doi:1 0.1 371/journal.pone.009781 8.g002 



Microarray database analysis for genome-wide 
methylation and expression assay 

The lUumina Infinium Methylation Chip data for breast, colon, 
liver, lung, and stomach cancer were obtained from the GEO 
database (http://www.ncbi.nlm.nih.gov/geo/). Microarray data 
for ten sets of normal tissue and its matching cancer tissue for each 
cancer type (fifty sets in total) were collected. The accession 



numbers of the adopted databases are indicated in Table S 1 . The 
lUumina Infinium Methylation Assay detects genome-wide meth- 
ylation covering 27,578 CpG sites in 14,495 genes. A methylation 
index (P) was outputted for each site, which is a continuous 
variable ranging between 0 (no methylation) and 1 (100% 
methylation), representing the ratio of the intensity of the 
methylated-probe signal to the total locus signal intensity. These 
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Figure 3. Pathways most strongly associated witKi the significantly altered genes in the five cancers. Top functional categories are given 
for breast cancer. The Ingenuity software assigns a p-value based on the likelihood of obtaining the observed number of category or pathway-related 
molecules in a given data set by chance alone. The threshold line denotes the p = 0.05 level. The line graph represents the ratio of affected genes to 
the total number of genes in a pathway. The top functional category for the other four cancers is added in Fig. S3. 
doi:1 0.1 371 /journal.pone.009781 8.g003 



values were then used to calculate a ratio of relative methylation 
between normal and tumor tissue, with higher values correspond- 
ing to greater levels of methylation in tumor tissue relative to 
normal tissue. 

In sUico expression analysis for selected genes was carried out on 
the Oncomine database platform (https://www.oncomine.org). 
To do this, datasets from the corresponding cancer type were 
collected for each gene, and expression was examined with a 
threshold satisfying p<0.01. 

Pathway analysis 

To identify pathways displaying methylation-specific altered 
expression patterns with potential roles in the five cancer types, 
functional categorization and pathway constructions were con- 
ducted using the Ingenuity Pathway Analysis (IPA) software tool 
produced by Ingenuity Systems as previously described [22]. A 
repertoire of genes from each cancer type, which was selected as 
differentially methylated genes was submitted to the IPA. The p- 
value for individual network was obtained by comparing the 
likelihood of obtaining the same number of transcripts or greater 
in a random gene set as are actually present in the input file (i.e., 
the set of genes differentially expressed in the normal tissue and 
cancer tissue), using Fischer's exact test based on the hypergeo- 
metric distribution. The highest confidence functional network 
was designated as the top network. 

Real-time RT-PCR 

Isolation of Total RNA from cell culture and reverse 
transcription were conducted as described previously [22]. 
Expression levels of selected genes were measured by real-time 
quantitative RT-PCR analysis. At least triplicate reactions were 
performed for each sample using a Kapa SYBR Fast qPCR Kit 



(Kapa Biosystems, Woburn, MA) with gene-specific primers 
(Table S2) on an ABI 7300 instrument (Applied Biosystems). 
RNA quantity was normalized to GAPDH content, and gene 
expression was quantified according to tire 2 ^'^^ method. 

Statistical analysis 

A hierarchical clustering dendrogram was generated by Gene 
Cluster 3.0 (http://bonsai.hgc.jp/~mdehoon/software/cluster/ 
software.htm) to determine the relationship between the ge- 
nome-wide methylation trend and cancer sample. In the analysis 
of microarray database, batch effects resulting from the use of 
different arrays at different tissues and time points were removed 
using ComBat [23]. Observations with adjusted p-values equal to 
or greater than 0.05 were removed and thus precluded from 
further analysis. Following adjustment, the remaining genes were 
defined as differentially methylated if [P (cancer) - fi (normal)] > 
0.3. A student's t-test was used to detect differences in the mean 
levels of methylation for selected genes, as well as for the 
expression level between the control and cancer cells. Statistical 
analyses were conducted using SPSS for Windows, release 17.0 
(SPSS Inc., Chicago, IL). A p<0.05 was considered to be 
statistically significant. 

Results 

In silico methylome analysis identifies tissue-specific 
nnethylation profile in breast, colon, liver, lung, and 
stomach cancer 

To identify and compare essential pathways that are deregu- 
lated by abnormal methylation of CpGs at the promoter in the 
major five cancers, breast, colon, liver, lung, and stomach, a 
repertoire of genome-wide methylation databases was first 



PLOS ONE I www.plosone.org 



4 



May 2014 | Volume 9 | Issue 5 | e97818 



Epigenetically Aberrant Regulatory Networks in Cancer 



GPX2 



CLDN15 



T0M1L1 



T3 

a> 



c 

03 
O 
I 

c 

(D 
T3 
03 

E 

CM 
D) 
O 



0.5 
0.0 
-0.5 
-1.0 
-1.5 
-2.0 
-2.5 
-3.0 
-3.5 
-4.0 
-4.5 I 




Breast 



Mixed Lobular 
and Ductal 
Breast Carcinoma 



03 



c 

0) 

u 

I 

c 

CO 
TD 
03 

E 

C\J 
D) 
O 




Breast 



Breast Carcinoma 



4.0 
3.5 
3.0 
2.5 
2.0 
1.5 
1.0 
0.5 



0.0 



T3 
0) 



C 
03 
O 
I 

c 
ro 

73 
03 

E 



Breast 



Mucinous 
Breast Carcinoma 



KL 



KISS1 



g 
"(5 



•D 

03 



C 
03 
O 
I 

c 
ro 

73 
03 



CM 
D3 
O 



6.0 

5.0 
4.0 
3.0 
2.0 
1.0 
0.0 
-1.0 
-2.0 



Breast 



Invasive Ductal 
Breast Carcinoma 



o 



2.0 
1.5 
1.0 
0.0 
-0.5 
-1.0 
-1.5 
-2.0 
-2.5 



Breast 



Mucinous 
Breast Carcinoma 



Figure 4. In silico expression profiles of selected genes with altered methylation in breast cancer. Five genes that have shown 
methylation change in breast cancer were selected, and their expression was examined in silico in the expression dataset and represented as a box 
plot. There appeared a coincidence that hypermethylated genes showed downregulation, or that hypomethylated genes showed upregulation in 
cancer. Data for the other four cancers is added in Fig. S4. 
doi:10.1371/joumal.pone.0097818.g004 



collected from the GEO database (http://www.iicbi.iilm.iiih.gov/ 
geo/) that harbors the lUumina methylation array data covering 
14,495 genes and 27,578 CpG sites. Datasets for ten pairs of 
cancer tissues and the corresponding normal tissues were collected 
for each cancer type except for breast cancer, in which case the 
data consisted of only one normal tissue sample and ten cancer 
tissue samples due to the lack of normal/ cancer parrs. 

To get insight into the relationship between the overall 
methylation change and each of the five cancer types, and further 
into the relative distance of the cancers, a cluster analysis was 
carried out. In the cluster, the normal and cancer tissues were 
generally separated in a bipartite mode, indicating that genome- 
wide methylation change is a key event distinguishing cancer cells 
from normal cells regardless of tissue types (Fig. SIA). Moreover, 
each cancer type was generally grouped together, thus explaining 
the tissue-specific methylomic profile in both normal and cancer 
cells (Fig. lA). A few cancer tissues were grouped nearby other 
tissues. For example, three colon cancer tissues were separated 
from the others and grouped together with lung cancers, showing 
hypermethylation in the VHL and PRDM15 (Fig. SIB). VHL 
encodes von Hippel-Lindau tumor suppressor and was shown to 
be hypermethylated in lung cancer [24]. The two subgroups of 
lung cancer tissues separated by the three colon cancers showed 
different cancer stages, with one subgroup including only stage I 
and with the other subgroup including also stage II and III. 



Interestingly, we found that breast cancer was closest with lung 
cancer whereas it was farthest from colon cancer. Next, we 
selected 90 key cancer-related genes and carried out the cluster 
analysis for cancer tissues to examine whether they could possibly 
represent the whole methylation status of cancer tissue. The results 
indicated a high ratio of overlapping between the two cluster 
analyses with the whole genome and the cancer-related genes 
(Fig. IB). 

Aberrantly methylated genes are involved In cancer- 
related pathways 

To construct the highest functional network in the five cancers, 
a collection of genes with differential changes of methylation are 
compared to the corresponding normal tissue, which was selected 
after filtering off the statistically insignificant genes (/)>0.05). In 
each cancer, the filtered genes were compared with another tissue, 
and genes appearing common in the compared tissues were 
selected. This comparison step was repeated using the ten normal- 
cancer pairs until less than 500 genes were filtered. At the end of 
filtration 3-4 pairs were included to construct the network. 

All of the sites that fit our significance criteria for differential 
methylation were examined for flinctional analysis using the IPA 
software tool. Of note, the "cancer" related network or bio- 
function appeared as the highest confidence among all the five 
cancers (Fig. 2, Fig. 3, Fig. S2, and Fig. S3). 
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Figure 5. RT-PCR analysis of selected genes that showed altered methylation in the five cancer types. The expression level of selected 
genes was examined using real-time RT-PCR in a normal and a cancer cell line for the five cancer types. In case of stomach, HEK293 was used as the 
normal cell line. There appeared a coincidence that hypermethylated genes showed downregulation, or that hypomethylated genes showed 
upregulation in cancer. Each sample was examined In triplicate, and the average expression levels are presented with standard errors. 
doi:1 0.1 371/journal.pone.009781 8.g005 



In breast cancer, ER-regnlated genes including FGFRl and 
IRSl, and ERK 1 -regulating genes including IRF6 and TOMILI, 
were aberrantly methylated. In addition, centered in the network 
is POU4F1, which interacts with ER and regulates INA, MSG, 
and CHRNA3. 

In lung cancer, two HOX genes, HOXA7 and HOXB4, were 
highly methylated. Meanwhile, FOSLl and THRB were hypo- 
methylated, taking centers in the network and appearing to 
function as the master regulator of a network of several gene 
transcripts potentially relevant to lung tumor development and 
progression. In rolon cancer, it is remarkable that AR is 
hypermethylated and acts as an interaction hub in the top 
network. Activation of AR in colon tumor is known to be a 
prosurvival signal, and it blocks migration [25]. PKA and ADC 
were also hypermethylated, which interacted with AR. In gastric 
cancer, VIM, which has been known to be a methylation marker, 
was also hypermethylated, and it centered in the network 
interaction with other proteins, such as PDGF complex and 
TWIST 1. In liver cancer, various kinds of GPCR, including 
ADRBl and DRDIPLC and their downstream element PLC, 
were abnormally methylated. 

Abnormally methylated genes are deregulated In cancer 

To examine whether the abnormally methylated genes induce 
deregulation of gene expression, in silico dataset analysis as well as 
real-time RT-PGR was conducted. Twenty-two genes, five from 
each cancer (two in case of liver), from the top hypermethylated 
genes were selected and in silico expression analysis based on 



Oncomine database platform was carried out. Among the 22 
genes 9 including KISSl [26], SERPINB5 [27], and GAMK2B 
[28] are previously established methylation markers in cancer or 
other diseases, meanwhile 13 genes including LILRB4, GLDN15 
are novel genes of which methylation status has not been known. 
As shown in Fig. 4 and Fig. S4, all the genes were dysregulated 
according to the methylation status in tumor cell lines with 
downregulation in case of hvpermethylation and vice versa. To 
further examine the relation between the expression and 
methylation, real-time RT-PGR analysis for five genes, one from 
each cancer, was carried out with RNA isolated from cultured cells 
composed of a pair of normal and cancer cell line for the five 
cancers. The result confirmed the fact that all the genes were up or 
downregulated according to their methylation status in the cancer 
ceU (Fig. 5). 

Discussion 

In this study, a key regulatory network was constructed from 
genes affected by promoter methylation in the highly incident five 
cancers. In addition, the relationship of the cancer types posed by 
the genome-wide methylation change- was elucidated. So far, 
cancer classification has been made based on the mRNA 
expression profiles [29-31]. Epigenetic changes have emerged as 
an alternative explanation to the signatures of cancer development 
in time courses as well as in different tissues. Recentiy, a 
quantitative atlas of histone modification signatures from 24 
human cancer cells revealed clusters of HEK293 and U20S; HFF, 
HeLa, and PC3; and MCF7, PANCl and MDAMB231. The data 
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showed a consistency between the mRNA levels and the histone 
modifying enzymes [32]. In another study, the cancer heteroge- 
neity across cancer types was addressed by the methylation 
variation in epigenetic domains, with tight clustering of methyl- 
ation levels in normal tissues, and marked stochastic variation in 
cancers [33]. These data suggest future efforts might instead be 
directed at defining the cancer epigenome as the departure from a 
narrowly defined normal profile. 

We found that breast cancer was closest with lung cancer and 
farthest from colon cancer. Colon cancer is notable in that most 
methylation alterations occur not in promoters, and also not in 
CpG islands, but in sequences up to 2 kb distant, the so called 
CpG island shores [34]. This unique epigenome characteristic 
might put colon cancer farthest from breast cancer. Breast and 
lung cancer have many common prognostic signatures, such as 
hypermethylation of BRGAl [35], SOX17 [36,37], and TLX3. 

It is interesting that the cluster that was constructed from the 90 
cancer-related genes showed a similar pattern at a certain level 
with that constructed from the whole genes, implying the cancer- 
related genes undergo a similar methylation change with that of 
the genome-wide change. It is therefore helpful to chase the key 
cancer-related genes when addressing the methylome status of a 
specific cancer. 

Notably, each cancer type represented a unique core pathway in 
the highest confidence network, implying that there appears to be 
a tissue-specific epigenetic cancer pathway. In the pathway, we 
found previously known epigenetic markers as well as novel 
markers that could be developed as epigenetic markers for cancer 
prognosis. In breast cancer, FLRT2 and SCTR are good 
candidates for novel markers; the expression of which is known 
to be down regulated [38,39]. However, no information for their 
epigenetic control is yet available. 

A limitation of our study is seen in the inability to assess 
methylation levels with respect to tumor grade, which thus 
precludes an examination of potential methylation patterns 
specific to cancer initiation versus those specific to tumor 
progression. In addition, detecting the frequency of methylation 
in tumor samples was limited due to the small number of databases 
for each cancer. Nonetheless, the degree of consistency between 
our findings and prior literature in the sampling of cancer-relevant 
genes gives bearing to the potential of using organ-specific 
alterations in methylation pattern as a prognostic tool in each 
cancer development and progression. 

In summary, we have presented a genome-wide methylation 
profile of five cancers, colon, liver, lung, gastric, and breast cancer, 
which we have compiled utilizing databases obtained from the 
most comprehensive methylation measurement technique cur- 
rently available. Our findings have provided information about the 
interrelationship among the five cancers with respect to methy- 
lome change. Moreover, we have found a key epigenetic 
regulatory pathway that is unique in each cancer. Further research 
into the mechanism leading to the common or dififerential 
mc-thylation in the cancers could give insights for estimating 
prognosis and determining treatment options. 



References 

1 . YokotaJ (2000) Tumor progression and metastasis. Carcinogenesis 2 1 : 497-503. 

2. Yun J, Frankenberger CA, Kuo WL, Boelens MC, Eves EM, et al. (2011) 
Signalling pathway for RKIP and Let-7 regulates and predicts metastatic breast 
cancer. EMBO J 30: 4500-4514, 



Supporting Information 

Figure SI Cluster analysis with gc-ncs of altered methylation for 
five cancers. Cluster analysis was carried out using a methylation 
index ((J-value) taken from lUumina Human Methylation 27 Bead 
Chip database registered at the GEO. A. Hierarchical clustering 
with the correlation distance was shown for 10 tumor tissues and 
10 normal tissues from each of the breast, colon, liver, lung, and 
stomach cancer. B. Significance analysis of microarray (SAM) for 
the 10 colon cancer tissues. Genes showing differential methylation 
among the 10 tissues were represented by the heatmap analysis. 
The color bar at top denotes the methylation level. 
(DOCX) 

Figure S2 Highest confidence network of genes displaying 
altered methylation levels in colon, liver, lung, and stomach 
cancer. In the network, hypermethylated genes in cancer are 
colored in red, whereas the hypomethylated genes are shown in 
green. The intensity of the color reflects the magnitude of 
methylation change. According to IPA, the network is relevant to: 
ceUular growth and proliferation, cancer, connective tissue 
disorders in colon cancer (A); post-translational modification, 
protein degradation, protein synthesis in liver cancer (B); cancer, 
dermatological diseases and conditions, cell death and survival in 
lung cancer (C); hereditary disorder, neurological disease, cellular 
assembly and organization in stomach cancer (D). 
(DOCX) 

Figure S3 Top functional pathways most strongly associated 
with the significantly altered genes in colon, liver, lung, and 
stomach cancer. Top functional categories are given for colon (A), 
liver (B), lung (C), and stomach cancer (D). The Ingenuity software 
assigns a P value based on the likelihood of obtaining the observed 
number of category or pathway-related molecules in a given data 
set by chance alone. The threshold line denotes the p = 0.05 level. 
The line graph represents the ratio of affected genes to the total 
number of genes in a pathway. 
(DOCX) 

Figure S4 In silico expression profiles of selected genes with 
altered methylation in colon, liver, lung, and stomach cancer. Five 
genes (two in case of liver) that have shown methylation change in 
each cancer were selected, and their expression was examined in 
silico in the expression dataset and represented as a box plot; colon 
(A), liver (B), lung (C), and stomach cancer (D). There appeared a 
coincidence that hypermethylated genes showed downregulation, 
or that hypomethylated genes showed upregulation in cancer. 
(DOCX) 

Table SI Genome-wide methylation databases analyzed in this 

study. 
(DOCX) 

Table S2 Sequences of primers employed in this study. 
(DOCX) 

Autlior Contributions 

Conceived and designed the experiments: SJK. Performed the experi- 
ments: BK SK GJ SEP SJK. Analyzed the data: BK SJK. Wrote the paper: 
SJK. 



3. Sabater L, Hoftberger R, Boronat A, Saiz A, DalmauJ, et al. (2013) Antibody 
repertoire in paraneoplastic cerebellar degeneration and small cell lung cancer. 
PLoS One 8: c60438. 

4. Easton DE, Bishop D T, Ford D, Crockl'ord GF (1993) Cjc-nctic linkage analysis 
in familial breast and ovarian cancer: results from 214 families. The Breast 
Cancer Linkage Consortium. Am J Hum Genet 52: 678-701. 



PLOS ONE I www.plosone.org 



7 



May 2014 I Volume 9 | Issue 5 | e97818 



Epigenetically Aberrant Regulatory Networks in Cancer 



5. Dick I'A, Rubin SM (2013) Molecular mechanisms underlying RB protein 
function. Nat Rev Mol CeU Biol 14: 297-306. 

6. Ai-MoundhriMS,Al-AnsariA,Al-MawaliK,Al-BahraniB(2013)BRCAl gene 
Molecular Alterations in Omani Breast Cancer Patients. Gulf J Oncolog 1: 45- 
51. 

7. MuUer I, Wischnewski F, Pantel K, Schwarzenbach H (2010) Promoter- and 
ccU-spccific cpigcnctic regulation of CD44, Cyclin D2, GLIPRl and PTEN by 
methyl-CpCi binding proteins and histone modifications. BMC Cancer 10: 297. 

8. Barckati Z, Radpour R, Kohler C. Zhang B, Toniolo P, et al. (2010) 
Methylation profile of 1T53 regulator) palh\va\ and mtDNA alterations in 
breast cancer patients lacking TP53 mutations. Hum Mol Genet 19: 2936—2946. 

9. OHveira AM, Ross JS, Fletcher JA (2005) Tumor suppressor genes in breast 
cancer: the gatekeepers and the caretakers. Am J CUn Pathol 124 Suppl: SI 6-28. 

10. Delgado-Cruzata L, Wu HC, liao Y, SanteUa RM, Terry MB (2013) Changes 
in DNA methylation by extent of breast cancer family history in unaffected 
women. Epigenctics 9. 

11. Xu Y, Diao L, Chen Y, Liu Y, Wang C, et al. (2013) Promoter melhyhilion of 
BRGAl in triple-negative breast cancer predicts sensiti\'ity to adjuvant 
chemotherapy. Ann Oncol 24: 1498-1505. 

12. Lee SJ, Chae YS, KimJG, Kim WW, JungJH, et al. (2013) AQP5 Expression 
Predicts Survival in Patients with Early Breast Cancer. Ann Surg Oncol. 

13. Carney WP, Bernhardt D, Jasani B (2013) Circulating HER2 Extracellular 
Domain: A Specific and Quantitative Biomarker of Prognostic Value in all 
Breast Cancer Patients? Biomark Cancer 5: 31-39. 

14. Wang YQ, Yan Q, Zhang JR, Li SD, Yang YX, et al. (2013) Ep^enetic 
inactivation of BRCJ A 1 through promoter hypermethylation in ovarian cancer 
progression. J Obstct Gynaecol Res 39: 549-554. 

15. Stefansson OA, Jonasson JG, Olafsdottir K, Hilmarsdottir H, Olafsdottir G, et 
al. (2011) CpG island hypermethylation of BRCAl and loss of pRb as co- 
occurring events in basal/ triple-negative breast cancer. Epigenetics 6: 638-649. 

16. Buckingham L, Penfield Faber L, Kim A, Liptay M, Barger C, et al. (2010) 
PTEN, RASSFl and DAPK site-specific hypermethylation and outcome in 
surgically treated stage I and II nonsmall cell lung cancer patients. Int J Cancer 
126: 1630-1639. 

17. Hino R, Uozaki H, Murakami N, Ushiku T, Shinozaki A, et al. (2009) 
Activation of DNA methyltransferase 1 by EBV latent membrane protein 2A 
leads to promoter hypermethylation of PTEN gene in gastric carcinoma. Cancer 
Res 69: 2766-2774. 

18. Montiel-Duarte C, Cordeu L, Agirre X, Roman-Gomez J, Jimenez- Velasco A, 
et al. (2008) Resistance to Imatinib Mesylate-induced apoptosis in acute 
lymphoblastic leukemia is associated with PTEN down-regulation due to 
promoter hypermethylation. Leuk Res 32: 709—716. 

19. Sebova K, Zmctakova 1, Bella V, Kajo K, Stankovicova I, et al. (2011) 
RASSFIA and CDHl hypermethylation as potential epimarkers in breast 
cancer. Cancer Biomark 10: 13-26. 

20. Seniski GG, Camargo AA, lerardi DF, Ramos EA, Grochoski M, et al. (2009) 
ADAM33 gene silencing by promoter hypermethylation as a molecular marker 
in breast invasive lobular carcinoma. BMC Cancer 9: 80. 

21. Liu Y, Song H, Pan J, Zhao J (2013) Comprehensive gene expression analysis 
reveals multiple signal pathways associated with prostate cancer. J Appl Genet. 

22. Kim SJ, KeUy WK, Fu A, Haines K, Hofiman A, et al. (201 1) (Genome-wide 
methylation analysis identifies involvement of TNF-alpha mediated cancer 
pathways in prostate cancer. Cancer Lett 302: 47—53. 



23. Johnson WE, Li C, Rabinovic A (2007) Adjusting batch effects in microarray 
expression data using empirical Bayes methods. Biostatistics 8: 118—127. 

24. Dmitriev AA, Kashuba VI, Haraldson K, Senchenko VN, Pavlova TV, et al. 
(20 1 2) Genetic and epigenetic analysis of non-small cell lung cancer with Notl- 
microarrays. Epigenetics 7: 502—513. 

25. Gu S, Papadopoulou N, Nasir O, Poller M, Alevizopoulos K, et al. (2011) 
Activation of membrane androgen receptors in colon cancer inhibits the 
prosurvival signals Akt/bad in vitro and in vivo and blocks migration via 
vincuhn/aetin signaling. Mol Med 17: 48-58. 

26. Moya P, Esteban S, Fernandez-Suarez A, Maestro M, Morcnte M, et al. (2013) 
KiSS- 1 methylation and protein expression patterns contribute to diagnostic and 
prognostic assessments in tissue specimens for colorectal cancer. Tumour Biol 
34: 471^79. 

27. Mardin WA, Petrov KO, Enns A, Senninger N, HaierJ, et al. (2010) SERPINB5 
and AKAP12 - expression and promoter methylation of metastasis suppressor 
genes in pancreatic ductal adenocarcinoma. BMC Cancer 10: 549. 

28. KimJH, Kim lAV, Kim SJ {201 1) Downregulation of ARFGEFl and CAMK2B 
by promoter hypermethylation in breast cancer cells. BMB Rep 44: 523-528. 

29. Patsialou A, Wang Y, Lin J, Whitney K, Goswami S, et al. (2012) Selective gene- 
expression profiling of migratory tumor cells in vivo predicts clinical outcome in 
breast cancer patients. Breast Cancer Res 14: R139. 

30. van 't Veer LJ, Dai H, van de Vijver MJ, He YD, Hart AA, et al. (2002) Gene 
expression profiling predicts clinical outcome of breast cancer. Nature 415: 530— 
536. 

31. Chou HL, Yao (TF, Su SL, Lee CY, Hu KY, et al. (2013) Gene expression 
profiling of breast cancer sur\'ivability by pooled cDNA microarray analysis 
using logistic regression, artificial neural networks and decision trees. BMC 
Bioinformatics 14: 100. 

32. Leroy G, Dimaggio PA, Chan EY, Zee BM, Blanco MA, et al. (2013) A 
quantitative atlas of histone modification signatures from humEin cancer cells. 
Epigenetics Chromatin 6: 20. 

33. Hansen KD, Timp W, Bravo HC, Sabunciyan S, Langmead B, et al. (2011) 
Increased methylation variation in epigenetic domains across cancer types. Nat 
Genet 43: 768-775. 

34. Irizarry RA, Ladd-Acosta C, Wen B, Wu Z, Montano C, el al. (2009) The 
human colon cancer methylome shows similar hypo- and hypermethylation at 
conserved tissue-specific CpG island shores. Nat Genet 41: 178—186. 

35. Paul I, Savage Kl, Blayney JK, Lamers E, Gately K, et al. (2011) PARP 
inhibition induces BAX/BAK- independent synthetic lethality of BRCAl- 
deficient non-small cell lung cancer. J Pathol 224: 564-574. 

36. Chimonidou M, Strati A, Tzitzira A, Sotiropoulou G, Malamos N, et al. (201 1) 
DNA methylation of tumor suppressor and metastasis suppressor genes in 
circulating tumor cells. Clin Chem 57: 1 169—1 177. 

37. Yin D, Jia Y, Yu Y, Brock MV, Herman JCi, et al. (2012) SOX17 methylation 
inhibits its antagonism of Wnt signaling pathway in lung cancer. Discov Med 14: 
33-40. 

38. Xu Yj Wei K, Kulyk W, Gong SG (2011) FLR F2 promotes ceUular proliferation 
and inhibits cell adhesion during chondrogenesis. J Cell Biochem 112: 3440- 
3448. 

39. Mayor R, Casadome L, Azuara D, Moreno V, Clark SJ, et al. (2009) Long-range 
epigenetic silencing at 2ql4,2 affects most human colorectal cancers and may 
have application as a non-invasive biomarker of disease. Br J Cancer 100: 1534— 
1539. 



PLCS ONE I www.plosone.org 



8 



May 2014 I Volume 9 | Issue 5 | e97818 



